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We consider multiple scalar fields coupled to gravity, with special attention given 
r> \ to two-field theories. First, the conditions necessary for these theories to meet solar 

C^ ■ system tests are given. Next, we investigate the cosmological evolution of the fields 

to see if these conditions can be met. Solutions are found in the dust era, as well as 
radiation and cosmological constant dominated epochs. The possibility of inflation in 
these theories is discussed. While power law growth of the scalar fields can yield the 
appropriate conditions to meet solar system constraints, these solutions are unstable. 
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1 Introduction 

Scalar-tensor (ST) theories are alternative models of gravity which provide a theoretics 
framework within which general relativity (GR) may be tested. Many important tests so fa 
have used the post-Newtonian approximation, and have relied upon corrections to dynamic 
in the solar system. Such tests force ST theories to limits where they greatly resemble GF 
For example, the current limit on Brans-Dicke theories|ll] requires uj > 5OO|0|, where uj ^ o 
recovers GR and uj ^ 1 would be a natural value to expect a priori. 

However, despite the strong resemblance to GR in our solar system, gravity may be rac 
ically different in other regimes. The gravitational field in our solar system is weak, severel 
limiting the parameter space of gravity tested. Theories which differ from GR in stron 
gravitational fields, but agree with solar system constraints, are needed to further test GF 
Recently, a class of ST theories with multiple scalar fields has been proposed[0. These tht 
ories can satisfy the solar system criteria to arbitrary accuracy, but still diverge from GR i 
other limits, for example in the strong field regime around binary pulsars. Thus, such theorie 
provide an important test of GR in a previously sparsely tested regime. 

Aside from their importance as generalizations of standard ST theories, multiple scala 
field theories have a second motivation. Such additional scalars coupled to gravity appear i 
Kaluza-Klein|^ and string theories [^] which seek to unify gravity with other forces. Tests c 
the correctness of GR thus provide constraints on such theories. Furthermore, ST theorie 
combined with Grand Unified Theories can provide adjustments to inflationary models of th 
early universe P] which allow the phase transition to complete in old inflation [^ and remov 
the fine-tuning of new and chaotic models^]. The inflationary universe not only solves severs 
longstanding cosmological problems, but also provides the only currently known source of see' 
density fluctuations which obey the magnitude and spectrum of the microwave background 
found by COBE0. While these "extended" versions of old inflation are severely constrainc' 
by the solar system tests|TO|, multiple ST models may offer successful conditions for inflatio 
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In this paper we consider a particular set of multiple scalar ST theories which are 
straightforward generalization of nonminimally coupled models. The one-field nonminimall 
coupled model is related to the JBD theory by a redefinition of the scalar field. Our startin 
action is thus different from those used previously §, where the action was taken after 
conformal transformation had been made so that the gravity sector appeared normal. W 
initially consider a model with an arbitrary number A^ of scalar fields, and derive the condition 
necessary to meet the solar system constraints. These conditions depend on the asymptoti 
values of the fields far from the solar system. Previous work used the ansatz that these value 
were constants which satisfy all the equations. However, the actual values of these field 
should be set by their cosmological evolution. 

In order to study their cosmological evolution, we first find a generic set of solutions whic 
exhibit power law behavior, and give exact solutions in the dust, radiation, and cosmologies 
constant dominated eras. Next, we specialize to a two-field model, which elucidates the mai 
features of the behavior. The stability of these solutions in time is then examined, and sma 
deviations are found to diverge from the exact power law behavior as the universe expands 
Hence, while solutions of multiple field ST theories are possible which meet the solar syster 
constraints and field equations for cosmological evolution, they are unstable. This findin 
calls into question the ability of these models to provide a true test of GR. We close wit 
comments on the possibilities of success of these theories and their potential role in infiatior 

2 N Field Models 

We begin by considering the action 

S = J d^x./^ {f{<j))R - GAB9"^<l>A,a<l>B,^ - 2f/(0) + IQttLJ) , (2.1 

where is an A^ component scalar field whose individual components will be denoted wit 
capital Roman letters, U is a. potential. Gab describes the kinetic coupling of the fields, L, 
is the Lagrangian for other matter and summation over A and B is implied. The scalar fiel 



and rescaling of the 0, Gab may be taken to be diagonal with entries of ±1 for each scala 
field. All other conventions are identical to those of Misner, Thorne and Wheeler ||11||. /(0) = 
G = I/k^, where G is Newton's constant, recovers GR with scalar fields, while /(</>) = C,4 
with (f) a singlet is the standard nonminimally coupled model. Making the field redefinitio: 
$ = ^(f)'^ for singlet places the action into standard Brans-Dicke form. 

Varying this action with respect to the metric gives the gravitational field equations 

ig^^u (-fR + GAB9''^(t)A,a<pB,p + 2^) - f,^, + Q^^Uf + f R^, - G AB(t>A,t.(t>B,u = S^T^u, (2.2 

with T^j, the energy-momentum tensor corresponding to L^. The trace of this equation is 

-fR+ GAB9''^(l>A,a<l>B,p + ^U + 3C/ = SttT. (2.3 

Varying with respect to gives the scalar field equations 

JaR + 2GABn(PB - 2Ua = 0, (2.4 

with a subscript A referring to the partial derivative d/d(f)A- In order to be viable, metri 

theories such as these must satisfy the solar system experimental constraints. Of these, w 

concentrate in the next section on the effects arising from first order space curvature (PPI 

parameter 7). 

2.1 solar system constraints 

To lowest order in the post-Newtonian approximation, the metric around a test body of mas 



M in the solar system is 11 



ds' = -(a-'^\ de +{l + 27^) [dx" + dy^ + dz% (2.5 

where 7 = 1 recovers the value in GR. In the JBD theory, 7 = (1 + a;)/(2 + a;) and curren 
tests from the Viking lander give u > 5000. 

We will use the (00) and trace equations to determine the expression for 7 in the theor 
of equation ( p.l|) . Note that 



To obtain this last approximation, we removed the first term by using the fact that V0 = 
0{M) and keeping only lowest order in M. Then, making use of equation (p7^), we find 

n/ = —G^lfAfBR + Ua/bG^I (2.7 

The matter source for the solar system bodies will be taken to be pressureless dust, wit 
delta function stress-energy terms. Making use of the representation of the delta function b 
V^(l/r), we therefore write 

1 „o fM\ ^ ^ 1 „o /M\ 



Again, to lowest order, the necessary curvature components are 



Roo = -^V'goo = -V (^— j , R = 2(1 - 2^)V' (^— j . (2.9 

The (00) equation then implies 

- (/ + fAG^lfsKl - 27) V^ (^) + /V^ (y) + f/ + UaG-A/b = 2V^ (^) , (2.1C 

where G is Newton's constant. If we take the case U = 0, which we will consider throughou 
the rest of this paper, and define 

G = /aGabIb, (2.11 

then we find 

G-2-f{f + C) = -2G. (2.12 

Similarly, the trace equation yields 

(2/ + 3C)(l-27) = -2G. (2.13 

Combining these two equations then gives 

, = I±^. (2.14 

This theory will be identical to GR for solar system tests if 7 = 1, which then implies 



As long as the coupling satisfies this last relationship, multiple ST theories can exactly repl 
cate GR in the solar system without being identical to GR. However, theories which satisf 



(|2.15| ) may still be quite different in other gravitational regimes. This behavior is in contras 
to single field ST gravity, where the solar system constraints force the entire theory to b 
indistinguishable from GR. 

Of course, 7 need not be exactly unity, but need only satisfy the current experimentc 
limits. Indeed, nature may dictate that 7 is not 1 but only some value close to 1, in whic 
case GR would fail. If this should be the case, both single and multiple ST theories coul 
meet the new constraint with a C that is a small but nonzero value, while higher order effect 
would cause a larger discrepancy with GR in other regimes. 

Whether the condition of equation ( p.l5| ) can be met depends on the values of the scala 



fields. Because Gab niay be diagonalized and rescaled, the vanishing of C implies either th 
/a are all in the solar system, or that some of the scalars have negative kinetic coupling 
Also, if C = 0, then equations ( |2.12| ) and ( |2.13D demand that / be equal to the gravitations 



constant, setting another condition. For a large number of models, including the one hen 
setting the background value of the fields to zero certainly meets these conditions. In man 
other models, such as that of 0, values constant in time will suffice. However, these values ar 
actually determined by the cosmological evolution of 0, to which we next turn our attentior 

2.2 cosmological evolution 

For cosmology, we use a spatially flat Robertson- Walker universe, with metric given by 

ds^ = -dt^ + a^{t) (dx^ + dy"^ + dz^) , (2.16 

where a{t) is the time dependent scale factor. From the (00) component of the gravitations 
field equations (|2.2|), we get 



m^f - \Gab'<1>a4>b + 3Hf-U = Svrp, (2.17 

where H = a/a is the Hubble parameter, p is the energy density, and an overdot indicate 
differentiation with resDect to time. The trace eauation eives 



where P is the pressure of matter, which is assumed to be a perfect fluid. Finally, the scala 
field equations, (|2.4|), yield 

n(f)A = -Ic^^lfsR + UbGjI = 0, (2.1C 

with the Ricci scalar given by 

R = 6H+ 12H'^ (2.2C 

One equation of the set ( p.l7| ) - ( p.l9| ) is extraneous due to the symmetries of the spacetim 
and the contracted Bianchi identities. 
dust era 

In the current dust dominated universe, with U = 0, the constraint for the solar syster 
(|2.15| ) must be satisfied. Contracting equation ( |2.19|) with /^ gives 

/aUcPa = -\fAGA\fBR = 0, (2.21 

a rather strong constraint on the behavior of the scalar fields in a dust universe. The limit 
imposed by the null summation ( p.21j ) will be more obvious when we consider a two field cas 



in the next section. 

To better understand this constraint, we consider particular choices for /(0) of the form 

/(0) = /, + FAB<pA(f)B, (2.22 

with fc and Fab constant. This model is a straightforward generalization of the standar 
one field nonminimal coupling. There are no added dimensional coupling constants. Only /, 
which corresponds to a modified gravitational constant, has dimensions. 

Because these equations are nonlinear, general solutions will be difficult, if not impossibk 
to obtain. We therefore look for power law solutions, by making the ansatz 



a(t) = ae(^-j , ^A = bA[^-j , (2.23 

where Cc, tc, 6a, P, and q are all constants. Although this power law assumption restricts th 



constant solution. For GR, we know that the cosmological expansion is power law, wit 
p being 2/3 and 1/2 in the dust and radiation regimes, respectively. Thus, physically, ou 
solutions allow easy comparison with the standard cosmological behavior. 



With the above simplifications, equation (|2.21| ) gives 

FABbAbB—j^iq"^ -q + ^pq) = 0. 

tc 



(2.24 



This equation is solved either by g = 0, g = 1 — 3p, or FAsbAbs = 0. The last condition cause 
/ to be constant, f = f^ = G. The gravitational constant does not evolve. The first conditio! 
g = 0, corresponds to constant fields, while the middle condition relates the evolution of th 
fields to that of the scale factor. 

We now examine the scalar field equations, ( p.l9| ), which imply 



bA{q^ - q + Spq) = G^lFBcbc{-Qp + 12p'). 



(2.25 



If q equals either or 1 — 3p, then the left hand side must be 0, giving p = or p = 1/2. Th 
first case corresponds to a static universe, in confiict with observation, and will be ignorec 
We thus retain only p = 1/2. 

Next consider the (00) equation, which for power law behavior becomes 



3p^ 
t2 



fc + F, 



t\2'? 



q^Gd ft\^'' , GpqFd /^^ ^g 



2t2 



tr 



t' 



l-Kpd f i_ 



-3p 



(2.2e 



(2.27 



with the definitions 

Fd = FAsbAbB, Gd = GABbAbB 

in the dust era. As noted from equation ( p.25|) , ifg = Oorg = l — 3p, then p = 1/2. Howevei 
plugging these values into equation ( p.26|) shows that the (00) equation cannot be solvec 
because the powers of time will not balance. Therefore, of the three possible solutions t 
equation ( p.24|) , only the last one, Fd = 0, has the possibihty of being consistent with bot 
an expanding universe and the other field equations. 
Now, equation ( p.25|) implies that 



Using the fact that both Fd and Gd are zero, the (00) equation siniphfies to 

3p2 Svrpd f^V^^ 



AC L (X.. \v(- 



(2.2C 



identical with the GR case. Thus, because of the cancellation necessary to meet the sola 
system constraint ( ^.151 ), the dust dominated universe with power law expansion resemble 
GR in all cosmological aspects. However, this does not mean that other cosmologies, such a 
a radiation dominated universe, necessarily resemble GR, as we shall see below. 
radiation era 

For a radiation dominated universe the energy density is given by 

pit) = pra\t), (2.3C 

with pr a constant, and P = p/3. Since the current universe is not radiation dominated, th 
solar system constraint, equation ( p.l5| ), need not apply. Because the stress energy tensor fc 
radiation is traceless, the trace equation. 



K \tr 



-6j9 + 12p2) q-^Gr ( t\^'i 6g(2g - 1)F^ / t\^'' ISpqFr f t ^'^ _ 



\ Liiy / fj \ fjiy / (j \, Liiy 



(2.31 
is most convenient to consider first. Here tr is a constant, and F^. and G^ correspond to th 
quantities ( |2.27| ) in the radiation era. 

Balancing the powers of time then requires that p = 1/2, exactly as in the standar 
GR case. With this value, the first term in equation ( p.31| ) vanishes, and the resulting tim 
independent part is solved by 



We then plug this result into the (00) equation, which is identical with equation (|2.26| ) excep 

with a right hand side 

npr ft\ -^P 
^"^ ' ^ . (2.33 



af \t. 



li q = 0, then the fields do not evolve, and 



This is similar to the Einstein GR case, except that the gravitational constant is shifted b 

When q = —3Fr/{Gr + ISF^), the (00) equation that results is solved by any one of thre 

conditions. These are 

Fr = ^ q = 

Gr = ^ q = -i (2.3S 

Gr = -6Fr => q = -\- 

This first case just reproduces the g = case above. The latter two conditions cause cance' 
lations such that 

-4, (2.36 



_6_ ^ »7r/v 2 
4ac2 at "' 



which is exactly the same as in GR. All of the above solutions also satisfy the scalar fiel 
equation. 

Much as with the dust case, these solutions strongly resemble GR, due to cancellations c 
the scalar field terms. If the scalars are constant during the radiation era, the gravitationt 
constant may be shifted. This may have observable consequences on nucleosynthesis, fc 
example. Deviations from exact Robertson- Walker behavior, caused by primordial black hole 
or density fluctuations, could also lead to regimes where the scalar fields play a more dynami 
role. 
cosmological constant era 

When the universe is dominated by a cosmological constant, the matter is given by 

P = Po = constant, P = —p. (2.37 

Such a situation arises, for example, in the inflationary universe[^], where po is the energ 
density of a scalar field either trapped in a false minimum or slowly rolling down a potentia 
The (00) equation again is similar to ( |2.26| ), except that the matter side is now given b 
Svrpo, and Fa and Gd are replaced with the appropriate Fq and Go in this regime. Examinin 
the powers of time, there are three, namely —2, 2g — 2, and 0. Unlike the previous matte 
conditions, where the right hand side had p dependence, these three powers cannot in genera 
be matchedlil. However, we may find approximate solutions to the general case by assumin 



that one of the two terms dominates gravity. 

If the fields are small, then the gravitational constant may dominate the field contributioi 
1/ K^ ^ Fab4>a4>b- Neglecting terms involving the then gives the standard GR equation 

3H^ = SttkVo, (2.3g 

which is solved by 



a(t) = aoexp 



-TTK^Poit - to) 



(2.3c 



This exponential growth corresponds to an inflationary universe. The scalar field equation 
give 

GH^Fab^b - Gab{4>b + SH<Pb) = 0. (2.4C 

Contracting with (pA and solving the resulting equation gives 



,, , , -SHGo ± \9H^Gl + 24:H^FoGo 
<Pa = bAc'^^'-'^K q = ^ ,^° — (2.41 

as a solution. Thus, for the small field approximation, the may grow, decay, or oscillate. 

If the are large compared to the gravitational constant, then a power law ansatz for th 
growth of the scale factor and fields gives for the (00) equation 

(sp'^Fo - ^Go + QpqFoj e'^-^ = Snpotl'. (2.42 

Equating powers of time demands g = 1, so the fields grow linearly, and remain dominan 
over the GR term. Combining the field equations with the time independent part of the (OC 
equation gives two possibilities: 

p = =^ Gq = — IGttpo^o 

The first case gives a static universe, while the second has the power of expansion dependin 
on the parameters of the theory. For p > 1, the second could give rise to extended typ 



are found. By using the (00) and cj) equations, these conditions can aU consistently be met if Gq = —4/ 



inflation]^ if the phase transition is first order, or soft inflation p] if the potential is slowl 
rolling of either the new or chaotic type. 

An exact solution can be found if the terms of power 2g — 2 in the (00) equation cance 
In this case, equation ( p.38D is again valid, and exponential solutions for the scale factor a{t 



and the are found. By using the (00) and equations, these conditions can all consistentl 
be met if Go = — 4Fo and q = H(—3 ± -\/3)/2. In fact, this exact solution is really a specie 
case of our first approximation, when the terms ignored exactly cancel. 

In the inflationary scenario, the potential of a scalar field acts as an effective cosmologies 
constant to drive the expansion. For Go > 2Fo, equation ( p.43| ) gives power law inflation. " 



the potential is of the new[|T2| or chaotic [|13| type, then this solution is a generalization of th 
soft inflationary scenario]^. In soft inflation, modifications of Einstein gravity can remov 
the fine-tuning of potential parameters found in regular gravity. The same situation can aris 
here, with the added advantage of avoiding the solar system constraints. 

More interesting is if the potential is of the old inflation|T^ type, giving rise to a first orde 



phase transition. Simple nonminimal coupling allows the phase transition to complete [0 
However, to avoid an overproduction of big bubbles of true phase, which would produce a 
excess of microwave background anisotropy, requires parameters which violate the solar syster 
constraint. Our theory can avoid this problem, and can do it with parameters which exactl 
reproduce GR in the solar system regime. Unfortunately, as shown in the next section, th 
models which reproduce GR are unstable in the dust era. 

3 Two Field Case 

Studying the A^ field case has the advantage not only of being general, but also of allowin 
compact notation. However, the interrelationship of parameters in these models can becom 
obscured. We next consider the simpler case of a two field model to better elucidate sue 
features, which include some stability problems. The action for two fields, in analogy wit 
equation (|2.1|) , is 



/(0,V') = ^ + /ll0' + /l2V^' + M^- (3.2 

A derivative cross term could also exist, but a proper rotation of the fields will eliminate thi 
term, so we set it to without loss of generality. 

The solar system constraint, equation ( |2.15D , now becomes 

{hlu + 4/i^r^)02 + 4/i3(/iir7 + h2Uj)H + [hlr] + Ahluj)^)'^ = 0, (3.3 

which may be rewritten as 

r]{h3iJ + 2hi(t)f + ooihcj) + 2/i2^)' = 0. (3.4 

Therefore, uj and rj must be opposite signs. We choose u to be positive. Furthermore, b 
rescaling the fields in the original action, these kinetic constants may be chosen to be cj = ] 
?7 = — 1, without any loss of generality. Thus, one important fact clearly shown in the tw 
field case is that multiple scalar tensor theories must have negative kinetic terms if they ar 
to agree with post-Newtonian tests of general relativity. 



From the relation (|3.4|) , one can also see the interplay between the kinetic terms uj an 
7], the coupling to gravity through the hi, and the values of the fields themselves. Agaii 
the asymptotic values of these fields are determined by cosmology, and specifically, by thei 
evolution in the present dust era. We now study that evolution to derive further constraint 
on the parameters of the theory. 

In a Robert son- Walker universe, the field equations become 

{2hi(f) + h3'^){3H + QH^)-{4) + 3H^) = 

{2h2il^ + h^){3H + 6H^) + {'4) + 3Hip) = 0. (3.5 

Again, we make the assumption of power law behavior, 

oW =..(-), ^ = 6(-). * = e(-). (3.6 

The field equations ( |3.5| ) then combine to give 



Substituting this relationship back into the solar system constraint ( |3.4| ) implies 

b = ±c, (3.g 

which also yields the relationship between the hi 

h3 = T{h + h2). (3.9 

These conditions also mandate that / = I/k^ in the dust era, as noted earlier for N fields 
The (00) and trace equations for the two field case become 

3H'f + 3Hf - ^(02 - ^2) = gvrp, (3.1C 

{6H + 12H^)f + 3/ + 9Hf + ^"^ - ^"^ = M{p - 3P), (3.11 

respectively. As in the A^ field case, power law solutions exist. In dust, the universe wi 
expand with power 2/3 just as in GR, and the same relation for the energy density, ( 2.29 
still holds. The power q at which the fields evolve is 



q=\\-l±\ll-^{h2-h,) |. (3.12 

The solutions for two fields in radiation and cosmological constant dominated universes als 
proceed in analogous fashion, and are easily obtainable from the more general A^ field case b 
substituting the appropriate expressions for F and G. 

4 Stability Analysis 

While the above solutions are exact, a crucial aspect of the system's behavior is whethe 
these solutions are stable. Only exceptionally fine-tuned solutions will start with exactly th 
initial conditions to fit the form of the above solutions. Further, any slight physical fiuctuatio 
can move the system away from the exact solutions, and the solar system constraints will b 
violated. Only solutions which return to the exact behavior when deviating slightly will likel 
satisfy the observable tests of these theories. We therefore examine the evolution of initiall 



We start by writing the dust solutions for the case b = c, hs = —{hi + /i2) as 

2 

a{t) = ao^-)' [l + ei(t)], 
m = 6(^)'[l + e2(t)], 

m = ^(^)'[i+^3(t)], 



(4.1 



where e, <^ 1 and q is given by equation ( p.l2| ) above. The scalar field equations then becom 



ihi 



4/l2 



g(g + l) 



£2 



-2(g + l)^-e2 



62 .. 2(/ii + /i2) eg 



ei 



+ [8{hi - h2) - 3g] ^ + 3{hi - h2)ii = 



q{q + l] 



|-2(g + l)f -e3 + ^^^^if^| + [8(/., - h2) - 3g] | + 3(/.,-/.2)ei = 



(4.2 
= 
(4.3 



while the trace and (00) equations give 



{12q^ + 6q+-){h-h2) + 2q^ 



"' "' + [Q{2q + l){h - h2) + 2g] ^^-^ 



t2 



+3(/ii-/i2)(e2-e3) + 



62^2 Vt, 



t\-2^ 



4ei 16ei 

7r + V + «^' 



(6g' + 7gV29)'' '' ■ ">-= 



< +(3,^ + 2,)fe-.3) + ^(i)""(f +«0=0. 



First consider the case of GR, with (f) = ip = Q. The trace equation then gives 

2 



This equation is solved by 



-ei + -e'l + 3ei = 0. 



^ ^1 A2 
^1- t +^2/3: 



(4.4 
(4.S 

(4.6 
(4.7 



with Ai and ^2 constant. The perturbation to the scale factor decays in time, and hence Gl 

is stable, as expected. 

Now, consider the perturbations of the scalar fields. Subtracting one scalar field equatio 

from the other gives 

2. 

£23 + -(g + 1)623 = 0, (4.8 

where £23 = £2 — £3 and equation ( |3.12| ) was used to remove the £23 term. This equation i 
solved by 



with B a constant. For the Q = —j case the perturbation grows logarithmically, and th 
scalar fields are unstable. When g 7^ — | there are two roots, one of which grows in time wit 
positive power. In general, the scalar fields will pick up a combination of these two roots, wit 
coefficients depending on initial conditions in the dust era. Thus, the perturbation e23 wi 
generally have a component which grows in time, and hence the scalar fields are unstable. 

That the solutions are unstable is not a surprise, for one of the fields has a negativ 
kinetic term. Such kinetic terms cause instability by allowing for infinitely many negativ 
energy states when the system is quantized. They also can lead to the classical instabilit 
observed here. In addition, our power law solutions include the case where the fields are zer 
far from the solar system, so that these values too will be unstable to perturbations as th 
universe evolves. Hence, while multiple ST theories can satisfy the solar system constraints 
the requisite values will be unstable, and hence are unlikely to actually occur. 

There are possible modifications to evade this instability. One of the motivations fc 
multiple ST theories is that they arise in string theories and other models which attempt t 
incorporate gravity with the other forces. The action used here may be just an approximatio 
of the full theory, whose other terms stabilize the evolution. Or, on a classical level, a mae 
or potential term for the scalars could freeze the at some minimum, thus preventing th 
runaway growth of small perturbations. Any such modification would have to be examine' 
in more detail to see if stability arises. However, given the problems associated with negativ 
kinetic terms, multiple ST theories that are consistent with GR in the weak-field limit seer 
to be unstable. 

5 Conclusion 

Multiple scalar ST theories offer the possibility of providing a new theoretical framewor 
in which to test GR. Because such theories can mimic GR in the solar system, where regula 
BD theories are severely limited, but yield substantially different results in other regimes, the 
are especially valuable. Our results complement previous workfOl, which used values constar 



work also note the need for negative kinetic terms, and point out the potential for problem 
with quantum and stellar stability. We have found still another problem. In theories like thij 
with negative kinetic energy terms, the cosmological solutions which meet the solar syster 
constraints are unstable as they evolve in time, and thus can not be maintained. 

Our time evolving solutions for the scalars are consistent with both solar system tests c 
GR and the cosmological field equations. If these solutions can be stabilized by some mean 
such as an explicit potential term or a more complete theory of gravity, then they will be 
useful probe of GR. However, without such stabilization, the only way to make these theorie 
agree with GR in the solar system is to choose constants that make the theory approximat 
GR not only in the post-Newtonian approximation, but also in other regimes. This situatio 
is in analogy with the single scalar JBD theory, and going to multiple scalars would have n 
advantage. 

Finally, we note that our solution in the cosmological constant dominated regime exhibit 
power law growth when the fields dominate, as seen in equation ( |2.43[ ). For power greate 
than 1, this yields inflation. The extended inflation scenario]^ used the JBD theory to creat 
power law inflation, in which case the first order phase driving inflation would eventuall 
complete, in contrast with the GR case. However, the JBD theory was inadequate, becaus 
values of the JBD constant which met the solar system constraint produced too many bi 
bubbles, and thus an inhomogeneous universe [|l^]. One might hope to use multiple scalar S' 



theories to give both power law inflation as well as meeting both the solar system constrair 
and the big bubble constraint. However, this scenario will not likely occur, as the dust value 
of the fields are unstable. Should a stabilizing mechanism be found, then further investigatio 
of the inflationary scenario in these models would be warranted. 
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